Nuclear Reaction Rates 


The rate for a given reaction depends on the number density of the 
reactants, Na and Nx, the rate these reactants encounter each other 
(i.e., their velocity) and on the cross section for the reaction 


number of reactions/nucleus X /unit time 


number of incident particles/cm7 /unit time 


Moreover, because of the electrostatic repulsion of the reactants 
(and because of atomic resonances) the cross section itself is a func- 
tion of particle velocity. Thus 


Tax = NaNxo(v)v 


Of course, the particles in a star will have a distribution of veloci- 
ties, hence in order to calculate the true reaction rate, we have to 
integrate over a velocity distribution, (v). In addition, if a and X 
are identical particles, the number of particle pairs will not be NZ, 
but only 1/2 this value. Thus the actual reaction rate will be 


Tax = (1+6ax)7* NaNx a a (v) v $(v) dv 
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where ax is the Kronecker delta, and Na is Avogadro’s number. 
For ø(v) we can substitute in the Maxwellian distribution, which, 
in the center of mass frame, is 
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where u is the reduced mass 
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Thus, the reaction rate is 
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The value (cv) itself depends on three factors: 
e The probability of overcoming the coulomb barrier 
e The probability of a quantum-mechanical interaction 
e Whether the reaction occurs near a nuclear resonance 


NON-RESONANT REACTIONS 


Before any nuclear reaction can occur, a substantial potential barrier 
must be overcome. For a typical nuclear size of 1.2 fm= 10^ !? cm, 
this barrier is equivalent to 


2 
vag PE 

R 
This temperature is orders of magnitude higher than that found 
in stars. Thus, for nuclear reactions to proceed, particles must 
tunnel through the potential barrier. From quantum mechanics, 
the probability of doing this is 
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Quantum mechanics also says that the cross section for two inter- 
acting particles is proportional to the de Broglie wavelength of the 
particles (since each particle sees the other as a smear over length 
A = h/p). Thus, 
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Thus, as long as there are no resonances, 
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where the astrophysical cross-section S(E) is a slowly varying func- 
tion that contains factors which are intrinsic to the individual nu- 
cleus. With this definition, the nuclear reaction rate, as a function 
of temperature (11.2) becomes 
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where b = 0.997,Z,A'/2 MeV!/? and A is the reduced atomic 
weight in a.m.u. 


We can evaluate (11.6) (or, at least, approximate its solution) using 
the method of steepest descent. First, note that the integrand is 
a sharply peaked exponential; by setting the derivative of the ex- 
ponent to zero, it is easy to show that the “Gamow peak” occurs 
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Next, let’s expand the argument of the exponential as a Taylor series 
about Eo, t.e., 


LA 


f (Eo) 
2! 


E b 


-IT pi f(Eo)+f (Eo)(E-Eo)+ 


f(E) = (E—Eg)^4- 


If we substitute for b using (11.7), the first term is simply 
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The second term then disappears (by definition of f (Ep)). That 
leaves the third term: 
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Now let A = 44/ EgkT'/3. The reaction rate given by (11.6) then is 
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Now let's look at the lower limit of the integral. For tempertures 
TX 10? K, one can show using (11.7) that Eg > A. This being 
the case, the lower end of the integral contributes very little to the 
result. We can therefore change the lower limit to negative infinity 
without much penalty, and make the integral analytic. Thus 
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which, after a bit of simple math yields 
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and 
(ov) = K (AZaZx) ' S(Eo) ?e* (11.11) 
where K = 4.5 x 1014 if S( Eo) is measured in erg-cm?, or 7.2 x 10719 


if it is given in the more common (at least for nuclear physics) units 


of keV-barns. 


The previous analysis assumes that S(E) is constant at S( Eo) over 
the relevant range of the Gamow peak. We can improve upon this 
approximation by expanding S( E) as a power series, and repeating 
the cross section calculation. If we substitute 


S(E) = S(Ey) + EN (E — Eo) 


in (11.6) and do (quite a bit) of math, we find that, to first order 
in 1/7, the result of (11.9) should be corrected by a factor of 
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What this says, is that instead of using S(Eo) in (11.9) and (11.11), 
we should define a variable Sp and use that instead. This new 
variable is the weighted average of the astrophysical cross section 
over the peak, which is given by 


So = S(Eo)G(T) = S(Eo) + > (5g) T 
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Note that, in this form of (11.9 - 11.11), the temperature depen- 
dence is entirely in the variable 7, which is proportional to T~!/°. 
Formally, 
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(where a = 42.49 (Z2Z2. A)'/ and S is in keV-barns), but we can 
see the temperature dependence better by taking the logarithmic 
derivative of (cv) from (11.11), i.e., 
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For typical central temperatures of stars (between 10° and 10’ K), 
vx 5. Nuclear reactions in stars are extremely sensitive to temper- 
ature! 
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Energy level diagram for nucleus of !! B. 


RESONANT REACTIONS 


The nucleus of an atom has energy levels similar to those of elec- 
trons. If the energy of an interacting particle coincides with the 
energy of one those levels, a resonance occurs, and the probabil- 
ity of interaction is boosted by several orders of magnitude. In 
this case, the off-resonance interactions (described above) become 
insignificant. 


To compute the reaction rate for resonances, consider first that the 
reaction described in (10.4.1) is actually 


db X SH SY AS (11.15) 


When the reaction occurs, particle a and X form a compound nu- 
cleus Z, which is in an excited state. Just like an atom, this excited 
state can be described via a series of quantum numbers. Also, just 
like an atom, an excited nucleus will, after a very brief period of 
time, decay to a lower energy state. This lower energy state may 
be Y + b. 


To estimate the strength of a nuclear resonance, let's first perform a 
semi-classical estimate of the maximum theoretical cross-section of 
nucleus X to particle a. First, let's say that the angular momentum 
quantum number of the Z* state is /. Classically, the cross section of 
X to a would be 7b?, where b is the impact parameter. In the realm 
of quantum mechanics, however, the impact parameter is given by 
the de Broglie wavelength of the particle pair, Å = h/p. Since 
angular momentum is a quantized variable, this implies that the 
impact parameter is constrained so that b = £A. Moreover, since 
angular momentum must be conserved, a collision which excites 
a nucleus to a state with angular momentum £ must taken place 
between the @ and £ + 1 orbital. So 


c =r {L+ 1) AY —n(A? = (20 + 1) 73 (11.16) 


This argument (in addition to being non-rigorous) ignores particle 
spins, which interact with orbital angular momentum vectorially. In 
practice, the theoretical limit to the cross section is 
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where J is the total angular momentum of the resonance state, and 
Ja and Jx are the total angular momenta of particles a and X. 


Next, let's consider the stability of state Z*. Through the uncer- 
tainty principle, there is a relation between a state's lifetime and its 
energy width, specifically, I = h/7. Thus, the probability of a state 
decaying via method i (as opposed to all other methods) is 
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Finally, consider the energy width of a given state. In exact analogy 
to electron orbitals, there is a finite energy width to a nuclear state; 
just like atoms, the width of the state is given by 
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where E, is the mean energy of the state. (Again, this is due to the 
uncertainty principle. The total energy of the state is indeterminate; 
we only have a probability estimate.) 


We can now estimate the nuclear resonance cross section. ‘This cross 
section involve three terms. The first is simply the energy width of 
the state. The probability of an interaction depends on the exact 
energies of the incoming particles: the closer to the mean energy 
of the resonance, the higher the probability of interaction, i.e., o œ 
f(E). Next, one must consider the probability of a actually decaying 
into state Z* as opposed to all other states. ‘This probability is 
simply [,/T. Finally, the probability of state Z* doing anything 
except decaying back into particle a (and therefore not producing 
any reaction) is (T — I'4)/T. Therefore 
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If we substitute for Å = h/p using p = v2u E, and if there is only 
one possible decay option (i.e., ly =I — D), then 
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(This is the Breit- Wigner single-level formula.) 


To estimate the rate of a resonance reaction in stars, c(E) can be 
substituted into the integral of (11.2). Over the width of a typical 
resonance (~ 1 eV), the Maxwellian distribution barely changes, 
hence we can adopt its value at the resonance, t.e., 
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For simplicity, we can again take the lower limit on the integral to 
negative infinity, t.e., 
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This gives us 
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where E. and (wy), = w are given in MeV, and the reduced 
mass A is given in atomic mass units. Again, there is an extremely 
steep temperature dependence, with 
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For typical energies of ~ 0.1 MeV, and temperatures of Tg ~ 50, 
pr 20. 


Å few points 


e In astrophysical situations, reactions sometimes change from 
one regime to another, depending on the temperature of the star, 
and the energy of the resonance. 


e Sometimes, it is very difficult to find all the resonances of a 
nucleus in the lab, and if a narrow resonance is missed, the quoted 
reaction rate will be substantially off. 


e Heavy elements have many more resonances than light elements; 
by the time you get into the third row of the periodic table, all the 
reactions proceed through (overlapping) resonances. 


e If particle a is a neutron, then there is no Coulomb barrier to 
overcome, and the cross section for capture is nearly independent 
of energy. Neutron reactions are only important for heavier ele- 
ments (AZ, 60) and are not energetically important. (Their rate is 
primarily determined by the number of free neutrons in the star.) 


e Some reactions (including the initial reaction in the proton- 
proton chain) involved the weak nuclear force. The temperature 
dependence of these reactions may be different from either of the 
forms presented above. (They depend on overlapping wave functions 
and the strength of the individual reaction.) 


Summary for Reactions 


When electron shielding is included, the nuclear reaction rate equa- 
tion becomes 
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aX = 1 a 12.2.1 
Tax = (14+ bax) A AS (ov) ( ) 
or, in terms of number density 

Tax = (1+ fax) Na Nx f (ov) (12.2.2) 


For non-resonant reactions, (cv) can be computed directly from 
S (Eo), (aE es Za, Zx, Aq, Ax, and T, while for resonant reac- 
tions, (ov) is derived from (wy), Ep, and T. 


Now let's simplify the terminology by first assigning 
Aax = (ov) (12.2.3) 


and then defining the lifetime of species X against reactions with 
particle a to be 

dNx Nx 
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The form of this equation is now identical with that associated with 
radioactivity or other decay phenomena. (Unlike radioactive de- 
cay, of course, 7, depends on the external environment.) From the 
equation above, 
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Note that the Kronecker delta disappears from the final lifetime 
calculation. According to (12.2.1), the reaction rate for identical 


particles has a factor of two in the denominator, but since each 
reaction destroys two particles, this factor cancels out. Note also, 
that by (12.2.4) and (12.2.5), the total lifetime of particle X from 


all reactions is j i 
-=>- (12.2.6) 
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